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We introduce a new family of one dimensional liquids, which we label as "Staggered-Liquids" , in the 
phase diagram of the ID Kondo-Heisenberg model. It encompasses three distinct spin gapped liquids 
and a new Luttinger liquid phase. A Staggered Liquid is characterized by gapless modes with "large 
fermi-sea" signature in the charge density wave (CDW) mode, and in that the superconducting order 
involves the near condensation of charge- 2e Cooper pairs with finite center of mass momentum. In 
particular, the conventional gapless 2kF CDW and fc = pairing modes are absent. In the process, 
we analytically derive the phase transition from an intermediate coupling spin gap phase to the 
strong coupling gapless LL phase of the Kondo-Heisenberg lattice model. 



I. INTRODUCTION 

Multicomponent one dimensional electronic systems, 
of which the one dimensional Kondo-Heisenberg (K-H) 
model is a particular example, exhibit new phases that 
were unanticipated in the earlier studies of the one di- 
mensional electron gas (IDEG). The one dimensional 
Kondo-Heisenberg model consists of a IDEG interacting 
with a Heisenberg spin-i chain via spin exchange inter- 
action. In the present paper we characterize the stable 
fixed points of this model. 

In particular limits of parameters, we obtain well con- 
trolled analytical solutions which enable us to enumer- 
ate and characterize the quantum numbers of all gap- 
less modes. Gapless modes are properties of the fixed 
point. Thus, our analysis lists the minimal set of stable 
fixed points in the global phase diagram of the Kondo- 
Heisenberg model. Surprisingly, we find that there is 
a common feature to all the fixed points in that CDW 
and pairing gapless modes are obtained with "unusual" 
wave numbers, and hence give the label "staggered liq- 
uids" to the family of fixed points. In particular, whereas 
previously catalogued liquid phases of the IDEG have a 
gapless charge 2e pairing mode at fc = 0, in a Staggered 
Liquid this mode appears at non-zero wave-vector. Put 
differently, in a Staggered Liquid phase the dominant su- 
perconducting order involves the near condensation of 
Cooper pairs with finite center of mass momentum. Simi- 
larly, there is no gapless CDW mode at wavenumber 2kF- 
Our results are summarized in tables-1,2 of section-II. 

In previous publications , we have already charac- 
terized two distinct spin gap phases (at weak coupling 
and at a "Toulouse point" value of parameters). In 
the present paper we add a new gapless Luttinger liquid 
(LL) phase (labeled "Staggered LL) which is obtained 
by going away from the Toulouse point towards stronger 
coupling. Interestingly, though its mathematical form is 
similar to commensurate-incommensurate transition, we 
find the phase transition is first order. To our knowl- 
edge, it is the first analytic derivation of the phase tran- 
sition from an intermediate coupling spin gap phase to 
the strong coupling gapless LL phase in the Kondo lattice 
model. Moreover, a third distinct spin gap phase is ob- 



tained by introducing weak attractive interactions to the 
staggered LL (and hence we name it "staggered BCS"). 

We emphasize that in this paper we limit ourselves to 
cataloguing the stable fixed points. We do not discuss 
the range of their basins of attraction and the validity of 
the solutions away from the quantitatively controlled lim- 
its of parameters. In other words, we defer to a future 
publication the discussion of how exactly to construct 
the phase diagram as a function of Kondo interaction 
strength for a given discrete Kondo-Heisenberg lattice 
model at given incommensurate filling and with spin ro- 
tation invariance. 

The paper is organized as follows: In section-II, we 
define the model and the order parameters. Our results 
are summarized in tables-1,2. (The derivation of these 
results is presented in the ensuing sections). In section- 
Ill, we review the weak coupling limit {Jk ^ Jh) spin 
gap fixed point solution j|] . In section-IV, we review the 
Toulouse limit {Jh <C Jk ~ Ep) spin gap fixed point 
solution, with some extended discussion of the Unitary 
transformation. In section-V, we derive the phase tran- 
sition to a gapless LL away from the Toulouse point to- 
wards stronger coupling {Jh <C -Ef <C Jk)- In section- 
VI we make some additional concluding remarks. In or- 
der to facilitate the reading of the paper, a discussion and 
bosonization representation of order parameters is given 
in an Appendix. 



II. KONDO-HEISENBERG MODEL AND ITS 
ZERO TEMPERATURE FIXED POINTS 

A. The model 

The Kondo-Heisenberg model (|l]) consists of two in- 
equivalent interacting chains; one is a one-dimensional 
electron gas (described by the Hamiltonian JJi^^G Q-j^ 
and the other an antiferromagnetic Heisenberg chain of 
localized spins 1/2, {fj}. The chains interact via a spin 
exchange interaction with an antiferromagnetic coupling 
constant Jk > 0. 
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Hk = 2Jk ^ fj ■ s{xj) 



(2) 
(3) 



where s{xj) = il!l^{xj)^^il!p{xj) is the electron gas spin 
density operator at position Xj of the local spin fj of 
the Heisenberg chain. We focus on the low energy and 
long-distance behavior of the electron's correlation func- 
tions by taking the continuum limit of the electron gas 
and linearizing the IDEG dispersion relation about the 
fermi points, zLkp, with corresponding right and left go- 
ing electron fields, Ra- and Lo-; 



Where a —1, |. 
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FIG. 1. Kondo- Heisenberg model 

The effective Fermi wave numbers (in the sense of the 
generalized Luttinger's theorem |^) for the IDEG and 
the spin chain are 2fci? and 2k^^''^ = Ti/h respectively 
(where h the distance between the local 

spins of the Heisenberg chain). It is assumed that the two 
systems are mutually incommensurate, and that 2fci? is 
incommensurate with any underlying ionic lattice. The 
continuum limit is taken for the IDEG while the Heisen- 
berg chain is initially left discrete (and remain so in some 
of the limit solutions derivations). Therefore, the totality 
of our analysis is rigorously valid in the limit 2kF 3> 7r/6 
(i.e., the number of electrons is much larger than the 
number of local spin-i moments per unit length). 

The IDEG spin currents are decomposed into forward 
and back-scattering parts; 

s{x) = ^ji{Xj)^ij0{Xj) 
= Js (x) -I- fis (x) 

where Js {x) ^Jsr (x) + Jsl (x), Jsb = ^Rt^ac'Ra' ; 
JsL = \L'^(^aa'Lcr' are the ferromagnetic {q — 0) spin 
currents of right- and left-moving electrons respectively, 
and 



(x) = e-^'^^^'riB. [x) -f e+'^'^'^^^ni (x) 

where nR = R'^^^^La' ] rlh — Lt'^2~^'^' ^''^ ^^"^ stag- 
gered magnetization {q = 2fci?) components of the IDEG. 

Due to the incommensurate electron filling, back- 
scattering interaction terms are irrelevant in the renor- 
malization group (RG) sense, and for our purposes may 



be dropped from the Hamiltonian which describe the 
fixed points. As a result, the spin and charge sectors 
decouple. 



H 



dx [He + Hspin] 



The charge sector is described by a Gaussian model ^ 



KeUl (x) + -^v, {d.,<^,f 



Since in this paper we are not interested in the effects of 
anomalous ID exponents {Kc ^ 1), we will set Kc — 1, 
unless otherwise explicitly stated. The subsequent anal- 
ysis and manipulations deal only with the spin sector 
fields. The Kondo exchange interaction reduces to 



H 



K ^ Jk^Tj ■ J (a 



(4) 



B. Order parameters and staggered correlations 

Study of the different stable phases of the Kondo- 
Heisenberg array begins with an analysis of the gapless 
excitations of the decoupled fixed point. From there, as 
usual, we sort the phases by determining which of these 
excitations become gapped, and which remain gapless in 
the presence of the (Kondo) couplings between the f DEG 
and the Heisenberg chain. Since our ultimate goal is to 
study the coupled system, we need also to consider the 
character of gapless excitations constructed of composites 
operators from the two subsystems. An extensive expo- 
sition of the order parameters is given in the appendix. 
Below, we note only the modes which are relevant for a 
spin gap system. 

In the spin gap phases, only spin-0 modes may be gap- 
less. Thus, we focus our investigation on singlet pairing 
modes (charge-2e, spin-0) and CDW modes (charge-0, 
spin-0). The corresponding usual IDEG order parame- 
ters are: 



1 



(^sp^^[R\l\+L\R^ 
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CDW 



(5) 
(6) 



Modes of composite nature are: A composite odd-parity 
singlet pairing 



SP 



r, 



(7) 



and a composite particle- hole mode, Oc-cdw, which will 
play a central role in the ensuing discussion 



c-CDW = niDEG ■ T. 



(8) 
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Upon evaluating the corresponding correlation func- 
tions Xii^yx') = {Oi{x)Oi{x')), we find gapless modes 
with power law correlations of the form. 

X^ {xj - X. ) = xo {xj - X. ) , (9) 

where xo ~ . The staggering factor (—1)"' in the cor- 
relation functions (||) is effectively modulating the usual 
power law correlations by the reciprocal lattice vector, 
of the spin chain {t, }. As a result, the gapless modes are 
found in unusual finite momentum values: The singlet 
pairing are with momentum ^ (and there is no fc = 
singlet pairing with charge 2e), and the gapless CDW 
modes are at momentum 

2k*p = 2kF + - (10) 



(and not at 2fci? as CDW in free IDEG). These are the 
defining characteristics of a "staggered liquid". 

Insight into the gapless modes properties is gained by 
considering the so called ry-pairing modes at momentum 

±2kF, 

VR = ^T^I (11) 
VL = L\Ll 

corresponding to right and left going singlet pairs. In 
a bosonization representation (Appendix-A) it is easy 
to see that the 77-pairing operators depend only on the 
IDEG charge sector fields. The charge sector is unaf- 
fected by the relevant part of the Kondo and Heisen- 
berg interactions in all the zero temperature fixed point. 
Therefore, the gapless 77-pairing modes always exist and 
carry momentum 2kp as in free IDEG. It is instructive 
to define operators 

^-i=(ryH + 7y,.) (12) 

^"'"^ ^ ^ {rjn - Vl) ■ 

(thought in themselves they do not carry a well defined 
momentum quantum number). We found that an inter- 
dependence of the gapless modes is established by the 
following operator identities; 



OsP^[OcDW,v'''n (13) 

O=[0cw,0 (14) 

= [Oa-cDW,v'''n (15) 

Oc-sp=[Oc-CDW,v"'"']- (16) 



Hence, the gapless wavenumbers of CDW and pairing op- 
erators are always connected by momentum 2fci?. 



C. Main results: Staggered liquids fixed points 

For the purpose of characterizing fixed points, the issue 
of counting gapless modes requires clarification. Since 
the 77-pairing modes are gapless in all cases were the 
charge sector is gapless (i.e., at all the fixed points of 
the Kondo-Heisenberg lattice model at incommensurate 
filling), the interdependence of modes (given in equations 
|l^ , |l6|) implies that formally only the CDW modes need 
to be counted, while the pairing modes Osp and Oc-sp 
are redundantly derived from combinations of CDW and 
ry-pairing operators. In spite of that, since common dis- 
cussions in the literature are done in terms of the usual 
pairing order parameters Osp and Oc-sp, we will list 
them in our tables below. 

We have found four distinct fixed points associated 
with different parameter values: 

1. Jk ^ Jh <C i?F ; a spin gap phase at weak- 
intermediate coupling 

2. Jk ^ Ep ^ Jh ; A spin gap "Toulouse point" 
phase ^ at intermediate coupling. 

3. Jk ^ Ep :s> Jh ; A gapless staggered Luttinger 
liquid phase at strong coupling. 

4. Jk ^ Ep Jh ; A spin gap staggered BCS phase 
at strong coupling, with additional weak attractive 
charge interactions. 

We comment that all the fixed point Hamiltonians (as- 
sociated with the above noted phases) which we derived 
are in fact spin rotation invariant, even though the bare 
interaction parameters were in some cases breaking spin 
rotation invariance. It is an example of the possibility 
that the ultimate zero temperature fixed point can posses 
higher symmetry than the original microscopic model. 
Yet, we emphasize again that such issues do not affect 
the validity of our analysis for cataloging the fixed points 
of the most general microscopic Kondo-Heisenberg model 
(with or without spin rotation invariance). 

In the tables below, we characterize the above noted 
fixed points in terms of momentum quantum number of 
their gapless CDW and pairing modes ("X" signifies that 
the particular mode is gapped). 

Spin gap phases fixed points: 





Tabic- 1: 


Spin gap phases 






OcDW 


Osp 


c-CDW 


Oc-SP 


weak 
coupling 


X 


X 


2kp 


b 


Toulouse 
point 




7r 
b 


2k p 


7T 

b 


Staggered 
BCS 


2kp 


IT 

b 


X 


X 



Gapless Luttinger liquid fixed points: 



3 



Tablc-2: Gaplcss Luttingcr liquid phase 



OcDW 


OsP 


Oc-CDW 


Oc-SP 


Staggered 
LL 


b 


2kp 


b 



Obviously, the gapless LL is characterized by having also 
gapless spin density wave (SDW) modes and triplet pair- 
ing modes. 

Whereas previously catalogued liquid phases of the 
IDEG have a gapless charge 2e pairing mode at fc = 0, in 
all of the above noted phases this mode appears at non- 
zero wave- vector. Similarly, there is no gapless CDW 
mode at wavenumber 2kF- 



III. WEAK COUPLING LIMIT (Ja < Jh) SPIN 
GAP FIXED POINT 



In the weak inter- chain coupling limit 
Jk ^ Jh, Ep. 

It is allowed to make further approximation by taking 
the continuum limit also for the Heisenberg spin chain 
(such approximation is not valid in the opposite limit 
Jk 3> Jh, which is discussed in section-IV). The local 
spin chain field is then also decomposed into the smooth 
(ferromagnetic) and staggered (antiferromagnctic) com- 
ponents; 



(17) 



(Note: we will consistently use the subscripts "t, s" to 
distinguish the spin chain fields from the IDEG fields). 

In order to distinguish contributions coming from var- 
ious interaction terms, we introduce distinct Kondo cou- 
pling coefficients for forward scattering (Jj) and mixed 
interactions (Jm)', 



Hk = Jf {Jr 



-JI) 



(18) 



-Jmi-iy .(J^ + Ji) 



The mixed interaction, of the ferromagnetic IDEG com- 
ponent with the staggered impurity component (i.e., the 
Jm (—1)"' n,- • (J_R + Jl) term) has naive scaling dimen- 
sion |, but the oscillating (— 1)"* factor, which acts as 
an effective extra derivative factor (dx), renders this 
term to be perturbatively irrelevant in the renormal- 
ization group sense with respect to the free Hamil- 
tonian, Hq. The forward current-current interaction, 
Jf {JtR + Jtl) ■ (Ji? + Jl), has scaling dimension 2 and 
is marginal relevant and leads to the opening of a spin 
gap [£|4|. (The term will prove to be essential for 
understanding the Toulouse limit solution in section-IV) . 
Therefore, at incommensurate filling in the weak coupling 
limit, the Kondo-Heisenberg Hamiltonian (pi) reduces to 



H^eak = n'+nl + J}jdx (J], + J£) • [3% + Ji) (19) 



Kill {x) + -V, (a,$,) 



K 



where Vr , Vg are the spin wave velocities of the Heisen- 
berg chain and IDEG respectively {vr = 7rJ-H_/2). For a 
detailed derivation of the gapless modes of model (|is|), 
we refer the reader to our paper The end results are 
quoted in the first line of table- 1. It is remarkable that 
only composite modes are gapless. 



IV. TOULOUSE LIMIT (Jh < Jk 
FIXED POINT 



Ef) SPIN GAP 



In the limit 



Jh <C Jk ^ Ep 



the intra-chain interaction, Jh, is small compared with 
the interchain interaction, Jk, and it is incorrect to take 
the continuum limit for the spin chain prior to accounting 
for the effect of the interaction Jk- For simplicity, since 
Jh ^ Jk, we will first take the limit Jh = 0. (We 
shall find that bringing back Jh <C Jk is an irrelevant 
perturbation, in the renormalization group sense, due to 
the spin gap of "Toulouse fixed point" phase). Thus, we 
model the local spins as initially independent, and leave 
the Kondo interaction in its discrete form. 



IDEG 




2Jk V r. 



j 



In this limit, effective interaction and coherence between 
the local spins will come about explicitly mediated by 
the itinerant IDEG (i.e., in a kind of RKKY which is not 
introduce by hand to the model as Jh)- 

Below, we review and discuss the "Toulouse point" 
derivation and results of For the purpose of cal- 
culating correlation functions, wc bosonize the IDEG 
fermionic fields 




-iy^le^(x)-<p„(x)] 



Where 0,{x) = J^^dx'U,{x'), and [n,{x'), (l)^{x)] = 
—i5{x' — x), a The anticommuting Klein fac- 

tors, {F„,Fa'} — 5(j,a', are needed for the proper anti- 
commutation of fermions with different spin. As com- 
monly done, we re-express the operators in terms of 
bosonic spin fields <ps{x) = -^['/'T ~ "^ili ^^^^ charge fields 
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(j>c{x) = "^[0T + "^i]' correspondingly defined mo- 
menta lis and He. 

The crucial new step which we introduced in Q| is to 
make a unitary transformation of the fields, 



U = cxp 



(21) 



j 

UV2^id,(j},) = V2^id,^s) - 27r^rf<5(xz - x) (22) 

I 

UT+e~'^'^'U^=T+. (23) 



U cos 



27r0,(j) U^={-iycoa V2^^4j) 



(24) 



In words, going across an impurity, the spin-phase v 27r0s 
is shifted by ±7r. i.e., transformed fields with opposite 
spins acquires a phase shift of Jo- = ±-|. It is reminis- 
cent of the unitarity limit scattering we expect from the 
low energy physics of the single impurity Kondo effect 
Therefore, we interpret the unitary transformation 
as going to a Kondo strong coupling basis. 

The resulting transformed Kondo lattice Hamiltonian 
is given in equation (|25|), 

^ (25) 

V 71 " 

3 

Jl 

na 



3 

.5]t;(-1)-'"cos[V2^0,(j)] 

3 

Ho = + H§ - (J. + AJ,) i i^J 



where, ^ ^ f dx \n^^ + {d^cp^f 



and 



(26) 



(27) 



In ( p5| ) we have introduced independent interaction coef- 
ficient Jz and J± for the Ising and spin-flip parts of the 
Kondo exchange interaction Jk- Hence, formally we are 
here examining a generalization of the Kondo-Heisenberg 
model (|lj) to non spin-rotation- invariant interactions. 

The transformed fields constitute the low energy spec- 
trum of Hq into which part of the interaction energy has 
been incorporated. The transformed fields are taking ad- 
vantage of the Ising part of the magnetic Kondo interac- 
tion at the cost of kinetic energy (due to twisting of the 
spin field 4>s {x)). These are originally high energy states 
of the bare free IDEG Hamiltonian, H^^^'^ . For the 
transformed fields to become the new low energy states 
due to interactions, it is clear that the Kondo interaction 
strength needs to be on the order of the IDEG band- 
width. To this effect, note the shift of the ground-state 
energy per impurity (irrespective of the existence of a 
spin gap) in equation (Eq); 



= -{2Jz-ttvf)—. 

It represents the absorption of a part of the Kondo in- 
teraction energy, — (equal to the gain from forming 
an Ising singlet) into the transformed free field Hamil- 
tonian (p^), at the cost of kinetic energy +^^- Hence, 
for strong enough interactions the transformed free fields 
have lower energy than the bare IDEG free fields, and 
therefore determine the low frequency correlations of var- 
ious order parameters. Thus, the Toulouse point solution 
is an outcome of finite "strong enough" interactions and 
cannot be reached by perturbative methods about the 
non-interacting basis. 

For a special value of the coupling constants, 



TTVp 



AJz = 



(29) 



(the " Toulouse point" ) , we are left with an exactly solv- 
able fixed point Hamiltonian 



H* = m 



J 



m + —Y. ^/(-l)' cos[\/2^^,(j)]. (30) 



The spin part of the fixed point Hamiltonian has a dis- 
crete sine-Gordon form, and therefore, a spin gap. The 
transformed spin fields which develope an expectation 
value are (r/(-l)J) ^ and (cos[\/2^0s(j)]) ^ 0. In 
calculating correlation functions, it is important to re- 
member the effect of the unitary transformations which 
lead to 



cos 



27r0s(2:)]cos[V2^</.s(a;')]) (-1) 



{i{x) is defined as the j-impurity site to the left of po- 
sition x). The bare impurity correlations (^'j^j''^ de- 
cays exponentially. However, the transformed impurity 
spin, fj = U'^tJU , exhibit staggered long-range order at 

T = 0, = const ■ (-l)(-'-j'); This non-local order 

parameter characterizes the coherent ground state. That 
is all the information needed for deducing the correlation 
functions of all order parameters, and thus determine the 
gapless modes as done in [Q and summarized in the sec- 
ond line of table- 1. 

We take this opportunity to elaborate on the signifi- 
cance of the fields transformation. The tt phase shift of 
the field V27r0s across an impurity site (|2^) , give rise to a 
staggered coefficient (—1)^ in the Hamiltonian (|30| ) since 



exp 



z27r^- 



1=1 



= (-If 



Note that the factor (—1)"' is effectively "counting" im- 
purities, and is obtained irrespective of the order of the 
bare {rf} themselves (imagine an Ising chain of {rf }, 



(28) there is a factor e"^^^ = — 1 per impurity). Indeed, the 
correlation function {TjTp) is short range. 
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It is interesting to trace back the relevant interac- 
tion in the Toulouse fixed point Hamiltonian in terms 
of the continuum limit of the Heisenberg spin chain 
(p^). Due to the additional (—1)"' phase factor, in 
the transformed basis, the relevant slowly varying in- 
teraction is now Jmiir • (Jfl + Jl)7 while the interaction 
Jf (—1)"' (Jt_r + JrL)-(Jfl + Jl) is now also rapidly oscil- 
lating and irrelevant!! Thus, Toulouse fixed point physics 
originates from the interaction J,„ that couples the con- 
duction electrons to the staggered component of the im- 
purity array, an interaction that is relevant only with 
respect to the transformed fixed point Hamiltonian, i7o, 
and was irrelevant in the untransformcd basis. This pos- 
sibility would be missed in the continuum limit if we had 
dropped the (— 1)"' • [3^ +3^) term at the out- 
set (as is usually done, e.g., in Q). The perturbative 
relevance of various interaction terms is changed after a 
transformation to the "proper" strong coupling basis of 
fields about which perturbative RG analysis is performed. 
The notation J„j is not accidental, and it is exactly the 
one which is responsible for the non-trivial fixed point of 
the two-impurity Kondo problem [pj] . 



V. STRONGER COUPLING {.Jh <. Ef <. Jk) 
STAGGERED LUTTINGER LIQUID FIXED 
POINT 

A. Phase transition away from the Toulouse limit 

In a previous paper [Q, we 
analyzed the commensurate-incommensurate (C-I) tran- 
sition in the charge sector at the Toulouse point, as a 
function of the filling factor, and found a phase transi- 
tion from an insulating phase (with both charge and spin 
gaps) to a conducting phase with only a spin gap. Here, 
we are interested only in the case of incommensurate fill- 
ing (for which there is no charge gap). In this section, we 
analyze the phase transitions in the spin sector by vary- 
ing the parameter values away from the Toulouse point. 
While maintaining the same incommensurate charge fill- 
ing factor. 

We investigate the phase transitions within the trans- 
formed fields Hamiltonian (|25|). The local stability of 
the Toulouse limit fixed point (J* = nvp) is ensured 
by the existence of a spin-gap. This is all that can be 
deduced from perturbative renormalization group calcu- 
lations. Thus, the phase transition can be established 
only via non-perturbative methods. Below, we determine 
analytically the finite parameter space region character- 
ized by the Toulouse fixed point solution, i.e., the zero 
temperature stability of the spin gap to finite deviations 
AJz = (Jz — TTVp) > away from the Toulouse line to- 
wards stronger coupling. We find a new electronic gapless 
phase beyond a finite distance from the Toulouse point. 

Treating the transformed impurity spins in self- 
consistent mean-field approximation, we replace them by 



their expectation values in the transformed Hamiltonian; 

U^HU ^ + H^+^Jz\I^J2 i^D c>rM^j) (31) 

j 



The spin sector of the Hamiltonian ( |31| ) has a form fa- 
miliar from the study of Commensurate-Incommensurate 
(C-I) transitions; 



1 



dx [Hi + id^(t>s - sf] 

+h / dxcos[/3(j3s (x)] 



(32) 



Where, P = V2tt, 



(S = (5osin(7) = AJ^-W-(Tf) 



(33) 



h = ho cos {-f) = Ji_ |(rj) 



c 



where c = -, and 

a ' 

((-l)^r/) 



- cos (7) 
i sin (7) 



(34) 



The general character of the phase transition in the 
Hamiltonian ( p^ ) is well known The system remains 
commensurate until |i5| exceeds a finite critical value S'^. 
Therefore, the Toulouse limit is proved to be stable over 
a finite range of parameter space, AJz 7^ 0. 

Yet, care should be taken to identify the exact nature 
of the transition and the character of the ensuing gapless 
phase. The groundstate of ( |3^ is determined by the field 
configuration that minimize the energ y. A s we shall see, 
the C-I transition in our Hamiltonian (|3l|) is unusual due 
to the fact that the parameters S and h are themselves not 
constants, but instead are dynamic fields which need to 
be determined self- consistently by the additional mean- 
field minimization condition on (tJ) and (tJ). There 
are, in principle, three possible ground state solutions 
for the Hamiltonian (p2|): 

• Phase-1: A uniform spin gap phase, identical to the 
Toulouse point solution, with no finite gradients of 
d^M^j) i.e., {cos[V2^Mj)]) 7^ 0, <(-l)^T/) ^ 0, 
{dx4>s) — 0, and hence, also (rj) = 0. 

• phase-2: A gapless incommensurate spin "soliton 
lattice" ground state with periodic step-like kinks 
in the 4>s(xj) field. In such a phase, {dx<ps) 7^ 0, but 
stiU (cos[V^(/)^(j)]) ^ 0, and both ((-1)Jt;) ^ 
and (r/) ^ 0. 
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• Phase-3: A free gapless SDW phase; {dx<j)s) 7^ 0, 
(cos[\/27r(/'s(i)]) — 0. In that phase (tJ) — and 

The name "soliton lattice" comes from the classical 
solution, which has long range periodic order. Quan- 
tum fluctuations turn the long range order into power-law 
correlations, and thus the quantum ground state should 
properly be termed a soliton liquid. Nevertheless, this 
does not change the qualitative distinctions (in terms of 
non-zero expectation values) between the various phases. 
For simplicity, I will discuss the phases in classical terms. 

To find the transition points between the phases we 
need to compare, for a given AJ^ = (J^ — ttvf) ^ 0, 
the ground state energies of the spin gap phase (phase- 
1) with the that of the gapless phases. The usual result 
for commensurate-incommensurate transition, where the 
parameters h and 5 are constant, is that the soliton lat- 
tice solution (phase-2) has lower energy than the SDW 
solution, and transition is second order. This is not the 
case here, due to the fact that the parameters h and S are 
themselves inter- dependent dynamic variables. Thus, we 
need to minimized the ground state energy with respect 
to both the soliton spacing, I, (as usually done) and also 
with respect to the mean-field parameter, 7. 

The resulting commensurate-incommensurate transi- 
tion in the transformed ID Kondo lattice Hamiltonian 
(^), is first order. The argument is the Following: Re- 
member that rj and tJ are non-commuting. Therefore, 
if there was a second order transition to the soliton lat- 
tice phase than at the transition point, 5 w ^\/Tio, both 
5 ~ (rj) ^ and h ^ |('^/)| 7^ are less by a finite 
amount from their respective maximum value, a-nd 
ho. The energy of the soliton lattice at the second order 
transition is equal to the energy of the commensurate 
phase with the same value of h, which is always less than 
the maximum energy of the commensurate phase- 1 (for 
which h = ho and i5 = 0). Thus, we establish that the 
commensurate-incommensurate transition is necessarily 
first order. 

But, what is the incommensurate phase? There is no 
closed expression for the soliton lattice energy away from 
the dilute limit (i.e., far from the putative second order 
transition). Yet, we can analyze the competition between 
phase-2 and phase-3 in the dense soliton lattice limit 
(when the distance between soliton centers is less than 
the single soliton width). In that limit, the commensu- 
rate energy contribution (due to ft. 7^ 0) is exponentially 
small, while the dx4>s term contribution is approximately 
linear in S < So- Thus (again in contrast with the usual 
case of constant coefficients S,h ^ 0), the dense soliton 
lattice energy is less favorable than the SDW phase-3, 
(in which h = and S — 60). 

The above argument leads to two possible scenarios: 
Either there is a sequence of two first order transitions 
(phase-1 — >phase-2 — >phase-3). Or, there is one first or- 
der transition (phasc-1 — >phase-3). We conjecture that 



the second possibility is the correct one, and hence the 
phase transition occurs at S'^"^'^"''^ = \/2ho, i.e., 

(AJ.),„*w = \/^ (35) 

V ^ O g (jj 

In conclusion, at a finite deviation i^Jz) critical from 
the Toulouse point towards strong coupling, there is a 
first-order commensurate-incommensurate transition in 
the spin field (ps, in conjunction with transformed im- 
purity spin flop transition from {(tJ) ^ 0, (tJ) = O} to 
{{Tj) = 0,(tj) 7^ 0}. The transition is from the spin 
gap phase-1, to the gapless SDW phase-3, with no soli- 
ton lattice region. 

B. Staggered Luttinger liquid - a new strong 
coupling phase 

The transition in the spin sector to the gapless SDW 
phase leads to a new state which we call "staggered- 
Luttinger-liquid" . The staggered-LL is expressed in 
terms of the transformed fermion fields, which have com- 
posite phase fields, 

R,{x) = UR„{x)U^ = i?,(a;)e+''"^-^-<-"^" (36) 
L^{x) = UL,{x)U^ = L,(x)e"''"^-^<-"^" (37) 

where a — zL^ is the electron spin, and is the impurity 
operator which can take values ±i, (so 2TTT^a — 

Same as for the Toulouse point, we calculate the corre- 
lation functions with respect to the spectrum of the trans- 
formed Hamiltonian. All the order parameters which in 
Bosonic form depend on the (f>s field, have staggered cor- 
relation functions (as defined in (^) irrespective of the 
impurity configuration {tj}. Since there is no spin gap, 
there are now also gapless spin density wave (SDW) and 
triplet pairing modes. In the bosonization representa- 
tion, both cos {^y/2TT(f>s) and sin (\/27r(/)s) have power law 
decay of correlations (with an added staggered factor 

The Tj order of the transformed impurity array re- 
quires further clarification. The inter-impurity interac- 
tions generated by integrating out the transformed IDEG 
degrees of freedom in the residual Kondo interaction, 
AJ/ J2j Tjdxj4>s, are long ranged (i.e., well beyond near- 
est neighbor interaction). Honner and Gulacsi sug- 
gest that the effective interaction is ferromagnetic, and 
thus at least conforms with strong coupling calculations 

C. staggered-BCS; A third spin gap phase? 

It is interesting to investigate what would be the form 
of a BCS pairing of a composite staggered-LL? i.e.. 
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we introduce a conventional weak attractive interaction, 
U < (e.g., due to phonons), to the IDEG Hamiltonian. 
The singlet pairing take the form 



OsP - 



1 



V2 



L^Rl + R^Li 



(38) 



1 

71 



[L^Rl+R^L^] = (-1)^'(") OsP 



The resulting pair correlation function is staggered (as 
defined in (@)), witli nodes at the Kondo impurity peri- 
odicity. It corresponds to a neg ative Josephson coupling 
across each Kondo impurity . We stress that the node 
in the pair correlation function due to negative Joseph- 
son coupling is a node in the pair center- of-mass motion. 
It should not be confused with a node in the relative pair 
state There is no gapless fc = pairing mode. 

On the other hand, all the composite modes are now 
incoherent! In order to see this, note that in the gapless 
staggered LL phase, the gapless composite pairing mode 



Cc-SP came from the component sin (v27r(/)s) rj (see 
eq,51 in the appendix). Due to the singlet pairing inter- 
action (cos (V27r0s)) 7^ 0, and thus the correlation func- 
tion (sin (•\/27r(/)s (a;)) sin (\/27r0s (x'))) is exponentially 
decaying. Moreover, as in the staggered LL, (tJ) ^ 

and thus the part ^i?|L|r^ — i?|L|T+^ of Oc-sp is also 

exponentially decaying. These results are summarized in 
line 3 of table- 1. 

Our analysis suggests a new possibility: An unconven- 
tional staggered-BCS pairing phase may arise out of a 
two step process, where the staggered-LL is a precur- 
sor to the staggered-BCS phase; First, at a temperature 
Tfif, set by the renormalized Kondo interaction, there is 
a cross-over to a staggered-LL phase, characterized by 
the unitarity limit phase shifts. Then, at a much lower 
temperature, Tc, a conventional ECS pairing mechanism 
(e.g., phonons) leads to the un-conventional finite mo- 
mentum BCS pairing state. The above demonstrates the 
importance of considering the cross-over effects, due to 
strong interactions, prior to the consideration of pairing 
mechanisms. 



VI. CONCLUDING REMARKS 

The main results of this paper are: (1) The identifi- 
cation of the staggered liquids family of fixed point, as 
summarized in tables-1,2. (2) Derivation of the phase 
transition from a spin gap phase at intermediate cou- 
pling to a gapless staggered LJL at strong coupling. (3) 



The commutation relations (13-161) which relate CDW 



and pairing modes. Below we make some additional com- 
ments on our results. 

At weak coupling, the Kondo-Heisenberg model con- 
sists of a free electron gas coupled to a spin density wave 
system. One would naturally expect a BCS mechanism 
leading to a state of A; = BCS pairing of conduction 



electrons mediated by spin waves of the Heisenberg chain. 
We find it quite surprising that such a state does not ma- 
terialize at any stable fixed point of the one dimensional 
problem. 

Previous numerical simulations in the strong coupling 
limit have found that the "dominant" gapless CDW 
mode (in a gapless strong coupling LL phase) is with 
"large Fermi sea" wavenumber 2kp. Pairing modes were 
never evaluated. Yet, from our commutation relations 
( |T3| , p^ it is imperative that the pairing correlations are 
staggered!! We comment that, following the analysis in 
this paper, it is important that numerical simulations 
will establish the existence of both Ocdw and Oc-cdw 
gapless CDW modes. The pairing modes then follow au- 
tomatically as we explained. 

The numerical simulations were performed in the ex- 
treme strong coupling limit on a particular lattice struc- 
ture for which our analytical methods are not rigorously 
valid. Therefore, it is important to establish whether 
the gapless strong coupling LL phase in the numerical 
simulations is identical to the one which we derived an- 
alytically by a phase transition from the Toulouse limit 
solution |2|. 
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VII. APPENDIX: DISCUSSION OF ORDER 
PARAMETERS 

The study of the different stable phases of the Kondo- 
Heisenberg array begins with an analysis of the gapless 
excitations of the decoupled fixed point. From there, as 
usual, we sort the phases by determining which of these 
excitations become gapped, and which remain gapless 
in the presence of the (Kondo) couplings between the 
IDEG and the Heisenberg chain. In order to facilitate 
the readability of the paper, we give below the explicit 
expressions of various order parameters: 



1. Density wave modes 

The low energy spin currents of the IDEG, s{x), can 
be decomposed into two parts; 



■{x) = X {x) + [n, {x) 6^2'^^- + H.c] 



where 



A,<T 



-i'X,o 



(39) 

(40) 
(41) 



are respectively the k — Q and the k — 2kp components 
of the SDW mode (charge-0, spin-1) of the IDEG. (The 
index X = R,L correspond to Right/Left going electron 
fields). 

The Heisenberg chain spin current, f)- , may be similarly 
decomposed into a fc = part, Jr, and a finite momentum 



= f part, {-lyflr (where ^ 
vector of the Heisenberg chain); 



is the reciprocal lattice 



Jt (xj) + {-ly rir (Xj) 



(42) 



For the density- wave excitations, we count only the 
number of finite momentum excitations. It follows by 
symmetry that, for finite momentum, if there is a gapless 
mode at momentum q than there is also a gapless mode 
with momentum —q. We count them as one mode. To 
summarize, the gapless spin-1 excitations of the IDEG 
and the Heisenberg spin chain, and the operator whose 
correlation function is most directly sensitive to it are 
listed in table-Al. 



Table — Al: Gapless SDW excitations 



operator 


wave 
number 




2kF 


fir 


h 



(43) 



The incommensurate IDEG has one CDW excitation 
(charge-0 spin-0) with momentum, 2kp, created by the 
operator 



OcDW 



\,<J 



(44) 



(45) 



The generalized Luttinger theorem [^ asserts that there 
must be a gapless CDW mode at 2kp = 2kF + j- It 
is realized by the existence of composite-CDW [^ order 
parameters which are formed by combining a spin-1 SDW 
of the IDEG with a spin-1 SDW of the Heisenberg chain 
into a composite singlet Oc~cdw, 



O 



c-CDW 



1 



n r+) + 



Js ■ Jr + Js ■ fir (-1)"' 

e'2fe^^n^ . X + h.c. 
[e'^'^^'^ris ■ rir + h.c.] {-if 



n'T^ (46) 



(47) 



To summarize, the non-interacting two-chain system of 
a Luttinger liquid and a Heisenberg spin chain has gap- 
less finite momentum CDW modes at three wave vectors 
(table- A2). 

Table— A2: Gapless CDW excitations 



operator 


wave 
number 


fig ■ fir 


2kp + ^ 


OcDW 


2kp 


' J T 


2kp 


J s ' rir 


h 



(48) 
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Note that the composite-CDW excitations at wave vec- 
tors J and 2kp + ^ are not independent, since they can 
be related through a muhiphcation by the IDEG Ocdw 
(which has wave vector ikp)- Thus, there are only three 
independent gapless CDW modes. 



2. Singlet pairing modes 

The charge-2e singlet pairing modes also require care- 
ful consideration. In addition to the usual A; = BCS 
even parity singlet pairing, 



O 



SP 



1 
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g+.v^e. cos 



(49) 



be counted. The concomitant "trivial" existence of the 
usual BCS pairing Osp and composite pairing Oc-sp 
mode should be implicitly understood. 

The operator Oc-sp is odd under spin inversion oper- 
ation (i?j ^ T"*", as expected for a 
singlet. Note that it's spin inversion parity is odd, even 
though the conduction electrons part is in triplet pair- 
ing. In that sense the order parameter is a composite 
singlet. The operator is clearly odd under space inver- 
sion operation, P (exchanging R and L). The composite 
singlet operator, Oc-sp, can be arrived at by taking the 
tiine derivative of the BCS singlet order parameter 

oc [Hk,Osp\ = Oc~sPi where Hk is the Kondo- 



dt 



Hciscnberg Hamiltonian (2^). Therefore, Oc-sp is odd 
under time reversal, or alternatively, has only odd — w 
dependence. The corresponding order parameter on a 
discrete lattice (e.g. on a zigzag ladder) is [?] 



we note also the existence of an 77-pairing mode at mo- 
mentum ±2kp, 



Or 



-SP 



2 ^ 



(50) 



corresponding to right and left going singlet pairs. 

As with the CDW modes, in addition to the singlet 
pairing modes of the IDEG, it is necessary to consider 
the composite singlet pairing, Oc-sPj (a product of a 
triplet pairing in the IDEG with a spin-1 mode of the 
Heisenberg chain) which turns out to be odd parity [0,0 , 



(51) 



R\L\r- 



r\l\ +r\l 



trt 



+2i sin (^V2Tr(j)s'^ 



(The factor (— l)"* is needed so that both odd and even j — 
sites will conform in the continuum limit representation). 

There is a qualitative difference between the commu- 
tation relation (|lj) and previous commutation relations 
of Oc-SP in the literature. As elaborated below, we gen- 
erated the composite singlet, Oc-sp, by a combination 
of 2k*p composite particle-hole mode {nji ■ t) and finite 
momentum k = 2kp singlet {rjp pairing): 



riR 

Oc-CDW 
Oc-SP 



02kF-SDW 

[vl, Oc~cdw 



j ''''' 



(54) 



Note that the above generation of the composite sin- 
glet pairing is different from the usual way in which it is 
generated PI using the ^ momentum composite particle- 
hole mode {Jji ■ t) and k — Q momentum singlet {Osp 
pairing): 



(Note: If we do not take the Klein factors carefully into 
account than the bosonized form of the singlet and triplet 
composite pairing is erroneously exchanged!). It can be 
decomposed into two momentum components: a uniform 
fc = composite singlet 



1 



and a k = J , i.e. a staggered composite singlet 



(52) 



1 



Ot"|r (x) = (R'aa^L^) ■ (-1)-' (53) 



The commutation relations ( |l3| - |lq ) relate to each gap- 
less CDW mode a corresponding gapless pairing mode. 
Therefore, formally, only the 7/-pairing modes need to 



O. 



Jr. 
Osp 

c-SP 



1 

V2 



trt 



R.L 



1^1 



= Osp,Jr-t 



(55) 



Jp • T is an interaction term in the Hamiltonian, which 
develops a non-zero expectation value, (^Jr- ^ 0, in 
the spin gap phase of the Kondo lattice Hamiltonian. 



Osp, Jr. ■ t 



= [OsP,H] is 



Thus, the relation Oc-sp = 
important for establishing the time reversal symmetry of 
Oc-SP as determined by the Hamiltonian. The ^Jr ■ 
operator is not one of the gapless modes. In contrast, 
Oc-CDW = np • r is a gapless mode in the spin gap 
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phase. Thus, our relation Oc-sp = [vl,Oc-cdw], is 
estabhshing the inter- dependence of gapless modes in the 
spin gap phase. 

The commutation relations (??) and (??) indicate 
that there must be some symmetry difference between 
the usual CDW (Ocdw) mode of the IDEG and 
the composite-CDW {Oc-cdw) mode, and that the 
composite-CDW cannot be used in combination with 
j^euen^Q construct a BCS singlet mode, A, (as can be 
done with the usual CDW). Clearly there is no difference 
in the global symmetry properties of the two CDW modes 
(this would have been a violation of the generalized Lut- 
tinger's theorem). The difference is in a relative internal 
symmetry of the two chain system; a tt relative spin ro- 
tation around the z — axis of the IDEG with respect to 
the Heisenberg spin chain. This effect is best seen from 
the Bosonized spin field dependence of the operators, 



(where subscripts "1" and "2" refer to the IDEG and 
the impurity spin chain respectively). A tt relative spin 
rotation around the z — axis is shifting V27r [6is — 92s] by 
TT and leaving (jiis unaffected.. Thus under this operation, 
which we label 7?.f (tt), 



From these transformation properties (|56| ) it is clear 
that the composite-CDW cannot be used in combina- 
tion with Ty'^^^to construct a BCS singlet mode, since 



under 7^f (tt) , A — > +A, while 7?™™Oc-cdw 



—r]^^'^"Oc-cDW- Hence, our final conclusion is that the 
Toulouse point phase and the weak coupling limit spin 
gap phase of the Kondo-Heisenberg model are distinct 
phases (as summarized in table-1). 

There is a simple physical interpretation for what is the 
distinction made by the TZ^~^'^^ (tt) symmetry operation: 
The composite-CDW {Oc^cdw) is actually constructed 
out of two spin-1 SDW modes which are coherently com- 
bined into a total spin singlet. Therefore, the mode is 
sensitive to the coherent relative phases of the spin fields 
between the IDEG and the Heisenberg chain, which is 
probed by the Uf (tt). In contrast, the "pure" CDW 
mode (Ocdw) of the IDEG is independent of any rela- 
tive state of the Heisenberg chain. 




A 

Ocdw 

Oc^CDW 



+rj 
+A 

+ OcDW 
— Oc~CDW ■ 



(56) 
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